Multiplication of solutions for linear 
overdetermined systems of partial differential 

equations 



A large family of linear, usually overdetermined, systems of partial 
differential equations that admit a multiplication of solutions, i.e, a bi- 
linear and commutative mapping on the solution space, is studied. This 
family of PDE's contains the Cauchy-Riemann equations and the cofactor 
pair systems, included as special cases. The multiplication provides a 
method for generating, in a pure algebraic way, large classes of non-trivial 
solutions that can be constructed by forming convergent power series of 
trivial solutions. 

1 Introduction 

In this paper we study a wide class of linear first order systems of partial differ- 
ential equations, that allow a bi-linear multiplication in the space of solutions. 
The simplest example is the Cauchy-Riemann equations. We know that two 
holomorphic functions, f = V + iV and g = W + iW, can be multiplied in order 
to produce a new holomorphic function fg = VW — VW + i(VW + VW). In 
terms of the Cauchy-Riemann equations 
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this multiplication can be expressed in the following way: two solutions (V, V) 
and {W, W) prescribe, in a bi-linear way, a new solution (VW — VW, VW + 



T^W 7 ). From the basic theory of holomorphic functions, we know that any solu- 
tion of the Cauchy-Ricmann equations can be expressed locally as a convergent 
power series of a simple solution with respect to the described multiplication. 

The Cauchy-Riemann equations provide the simplest example of a system 
of PDE's that has a multiplication on its solution set, but there are more so- 
phisticated examples. One such example is the multiplication of cofactor pair 
systems, discovered by Lundmark in [5]. 

A cofactor pair system (or bi-cofactor system) is a dynamical system q h + 
T^cfcf = F h , h = 1, . . . , ii, on a (pseudo-) Riemannian manifold, such that the 
force F — F(q) has two different cofactor formulations F(q) — (cof J)~ 1 W = 
(cof J) -1 VV", where J and J are independent special conformal Killing tensors 
of type (1,1), V and V are smooth real- valued functions, cof J = (det J) J -1 , 
and V is the gradient ((Wy = g^djV). Cofactor pair systems have several 
desirable properties, in general they are completely integrable, they admit a 
bi-Hamiltonian formulation, and they are equivalent (or correspondent) to sep- 
arable Lagrangian systems [H HI El 13 [9] . 

A cofactor pair system is characterized by a pair of functions V and V , and 
a pair of special conformal Killing tensors J and J, that satisfy the relation 

(cof J)~ 1 W = (cof J)~ l VV. (1) 

For fixed special conformal Killing tensors J and J, the equation |T]) constitutes 
a system of first order linear PDE's for two functions V and V. In [5], Lundmark 
found that the equation ((T|) allows a multiplication of solutions. When n = 2 
the multiplication formula is given by 

(v, v") * (w, w) = (yw - det (j~ 1 j)vw, vw + vw - tr( ~r x J)V wj , 

where (V, V) and (W, W) are solutions of ([1]). We see that when det (J -1 J) and 
tr(J _1 J) are not both constant, we can choose trivial (constant) solutions (V, V) 
and (W, W) of ([T]) and obtain non-trivial solutions through the multiplication. 
When n > 2 a multiplication also exists, but one has to consider the related 
parameter-dependent system 

(cof (J + /zJ)) 1 W M = (cof J)~ 1 VV, 

which can also be written as 

(j~ x J + fj,I)VV^ = det (J~ 1 J + fJ,I)VV, 

where V^ is polynomial in the real parameter /x (note that throughout this pa- 
per, we use the notation V^ rather than V(fi) to indicate dependence on the 
parameter /i). 
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The most interesting property of this multiplication is that it provides a tool 
for producing new cofactor pair systems from known ones. Especially, infinite 
families of separable potentials can be constructed. For example, the Jacobi, 
Neumann, and parabolic families of separable potentials are all constructed in 
[TO] through a recursive process that is a special case of the multiplication of 
cofactor pair systems. 

Remark 1. In JSj, equations of the form ([1]), considered on a real or complex 
vector space where J and J are constant matrices, are studied, and the general 
analytic solution is described . 

In order to gain better understanding of this multiplication, systems of the 
form 

{X + fj,J)VV^ = det(X + /i/)W, (2) 

defined on a general (pseudo-) Riemannian manifold, were studied in 0], without 
referring to any underlying dynamical system. By analyzing the corresponding 
equations at each degree of fj, in the equation it becomes obvious that the 
equation is satisfied if and only if the degree of is n and the left hand side 
(X + ill) VV^ can be written as a product of the scalar det (X + fil) and some 
1— form which is constant in /i. We can therefore rewrite the equation @ as 

(X + /i/)VV^ = (mod det (X + /i/)) . (3) 

It turned out that the system (J2]) allows for a multiplication of solutions, similar 
to the one existing for cofactor pair systems, if and only if the tensor X satisfies 
the equation 

(X + /i/)Vdet (X + fil) = det (X + fiI)Vtr(X + fil). (4) 

Several classes of solutions of Q where discovered, and it became apparent 
that systems of the form (TT]) and the Cauchy-Ricmann equations only consti- 
tute special cases of a much larger family of systems of PDE's that admit a 
multiplicative structure on the solution space. 

It was also remarked in [1], that by considering more general systems than 
one finds other new classes of systems that allow multiplication. In this 
paper, we will examine that subject. The linear systems of PDE's that we 
consider are in general impossible to solve, but the multiplication provides a 
non-trivial superposition principle (on top of the ordinary linear superposition) 
that, given to solutions, prescribes a new solution in a bi- linear and pure alge- 
braic way. With this superposition principle, large classes of new solutions can 
be generated from known solutions. In particular, we can construct non-trivial 
solutions by forming convergent power series of a simple solution. The question 
then arises for which systems of linear PDE's these power series constitute all 
solutions, like in the case of the Cauchy-Riemann equations where all holomor- 
phic functions admit a power series representation. Besides providing us with 
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more systems of PDE's that admit a multiplicative structure on the solution 
set, the generalization helps us to better understand the multiplication for the 
systems already known (in particular the puzzling multiplication of cofactor pair 
systems). 

This paper is organized as follows. In section [3 we formulate an abstract 
framework for characterizing the class of systems of PDE's that admit multipli- 
cation. We define the *— operator and give a characterization of those systems 
that admit *— multiplication on the set of solutions. The multiplication pro- 
vides a method for generating, in a pure algebraic and non-trivial way, new 
solutions from known solutions. A second formulation of the systems, using re- 
lated matrices, is introduced. In this matrix notation the Euclidean algorithm 
for polynomial division, which is closely related to the *— multiplication, can be 
encoded in an explicit polynomial of matrices. Some algebraic properties of the 
multiplication are also mentioned in this section. In section[3]we investigate the 
explicit forms of systems that admit multiplication. The study splits into differ- 
ent cases depending on relations among certain discrete parameters that appear 
in the studied class of systems. Especially, some typical (generic) systems with 
multiplication are derived and examined. The most interesting property of the 
multiplication of solutions is that we can construct large classes of solutions 
by forming power series, with respect to the *— multiplication, of trivial solu- 
tions. Section|3]is devoted to study such power series solutions. The problem of 
constructing systems with *— multiplication is in general quite complicated. In 
section [5] several methods for constructing systems with *— multiplication are 
described. The last section |6] contains concluding remarks and natural questions 
raised by the study presented in this paper. 

2 Multiplication of solutions for linear systems 
of PDE's 

Let Q be a n— dimensional differentiable real manifold. Consider equations of 
the form 

A^dV^O (modZ^), (5) 

where is a (1,1)— tensor depending polynomially on the real parameter /i, 
and Z^, are real- valued C 1 functions on Q that also depend polynomially 
on fi. The expression A^dV^ is a 1— form which components are polynomial 
in fj,, and the unknown function is a solution of the equation ([5]) if these 
components are all divisible (when considered as polynomials in fj,) by the fixed 
function Z^. Let Z^ = Zq + [iZ\ + • • • + /i m_1 Z m _i + /i m , be a polynomial of 
degree m, then there is no restriction to assume that A^ = Aq+/j,Ai + - ■ - + iJ, k Ak 
has degree at most m — 1 (otherwise we can reduce it modulo Z^). In order to 
simplify the description of systems admitting a multiplication, we also assume 
that Vfj, = Vo + /iVi + • • • + /i m_1 y m _i has degree m — 1. 
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Remark 2. The system (J5]) should be compared with system We see that 
([5]) generalizes ([3]) in several ways. First of all, no metric is specified on the 
manifold Q corresponding to the system §5§ , and we consider a system of equa- 
tions for 1— forms using the exterior differential operator d, rather than a system 
of equations for vector fields expressed with the gradient operator V. Moreover, 
in §3§ the polynomial X + fil is of degree one and contains only one arbitrary 
tensor X , while in ([5]) we consider a more general polynomial A^, and instead 
of det (X + ill) , we consider an arbitrary polynomial that does not have to 
be related to A^. 

Since the highest order coefficient of Z^ is a unit, for each polynomial P M with 
coefficients in the commutative ring of real-valued functions on Q, there exists 
unique polynomials and R^ such that P M = Q^Z^ + R^, where degi? M < m. 
Thus, for each function V^, there exists unique 1— forms Bq, ■ ■ ■ , B m -i such that 

ApdVii = B + /iBi + • • • + /x m_1 -B TO _i (modZ^), (6) 

and the equation (JSJ) can be written as Bq = Bi = • • • = B m _i = 0. Thus, in lo- 
cal coordinates q 1 , . . . , q n , the system ([5]) constitutes a, usually overdetermined, 
system of nm first order linear partial differential equations for m dependent 
variables. We will see that there exist non-trivial systems ([5]) admitting a mul- 
tiplicative structure on its solution set. 

Define a bilinear operation *, on the set of all real- valued functions on Q that 
are polynomial in /1, by letting * be the residue of the ordinary product 
VpWfj. modulo Zp. In other words, * is the unique polynomial of degree 
less than m that can be written as V^W^ — Q^Z^, for some polynomial Q^. For 
certain choices of A^ and Z^, the *— multiplication maps solutions of ([5]) to new 
solutions: 

Theorem 1 (*— Multiplication). Let S denote the solution set o/((5J. Then * 
is a bilinear operation on S if and only if 

Af^dZft = (modZ^), (7) 

i.e., if and only if Z^ — fj, n G S. 

Proof. Given two solutions , £ S, let Q M be the polynomial such that the 
product VfjWfj, can be written as V^W^ = Q^Z^ + * W M . Then, we have 

= -Q^A^dZ^ {mod Z^). 

Thus, we see that A^dZn = is a sufficient condition for the existence of the 
bi-linear operation * on S. To see that it is also a necessary condition, consider 
the trivial solutions — fi and — fi" 1 ^ 1 . For this choice of solutions, the 
polynomial becomes a non-zero constant, which forces the relation ([7]) to 
be satisfied in order for * W M to be a solution. □ 
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The following algebraic properties of * are immediate consequences of the 
corresponding properties of multiplication in general quotient rings of polyno- 
mials: 

Corollary 2. The solution set S together with the scalar multiplication, addition 
(defined in the obvious way) and multiplication * is an algebra over R, where 
the * -multiplication is associative and commutative. 

To calculate the *— product * W^, we form the ordinary product 

V^W M = VqW + (V W 1 + Wo) (J.+ --- + 

{V m - 2 W m -i + V m -iW m - 2 ) v 2m ~ 3 + V m _ 1 W m _ 1 /i 2m ~ 2 , 

and replace fx m , . . . , /i 2m ~ 2 with their residues modulo : 

fj, m = —Zq — Zifi — ■ ■ ■ — Z m -i(j, m 1 

^rn+l = ZoZ m -i + (Z\Z m -\ — Zq) /! + ■■■ + (Z m -\Z rn -\ — Z m _ 2 ) //" 1 

In general it is hard to find non-trivial solutions of the system (J5]), but having 
the *— operator we can generate (in a pure algebraic way) an infinite family 
of non-trivial solutions by starting with trivial solutions. For example, we can 
construct non-trivial solutions by forming polynomials, or convergent power 
series, of the trivial solution fi 6 S : 

Vji = y j a r ^l, where /4 := /i* •■• (8) 

r factors 

with a r being real constants. We note that, since all trivial solutions can be 
expressed as polynomials in [i, every solution that is a sum of products of trivial 
solutions has the form 1(5)) again. As long as we consider domains in the manifold 
Q where Zq ^ 0, we can also allow negative powers /i~ a in ([5]), for every natural 
number a, by defining fi~ a ;= (/i _1 )° and 

M _1 := ~ ( Z i + ■ ■ ■ + Zm-iH m - 2 + A 1 ™ 1 ) , 

so that /i * = 1. 

The following two examples illustrate how the Cauchy-Riemann equations 
and the cofactor pair systems can be considered as special cases of systems of 
the form ([5]) that admit a *— multiplication of the kind described in theorem [1] 

Example 1 (Cauchy-Riemann equations). Let Q be the 2— dimensional Eu- 
clidean space with Cartesian coordinates (x,y), Z^ = 1 + /i 2 and 



-1 

1 
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where I is the identity matrix. Then the system ([5|) reduces to the Cauchy- 
Riemann equations, and the ordinary multiplication of holomorphic junctions 
follows from the *— multiplication: 

{v, v) * (w, w) = (vw - VW, VW + VW^j . 

Every solution (V, V) that satisfies the Cauchy-Riemann equations in the origin, 
can be expressed in a neighborhood of the origin as a power series of the simple 
solution [x, y) : 

oc 

where (x+yfi)* again denotes the r 'th power with respect to the * — multiplication. 

Example 2 (Multiplication of cofactor pair systems). Let m = n and suppose 
that the tensor is linear in pi and has the identity mapping as the highest 
order coefficient, i.e. = X + fil. The system (J5J) can for this special case be 
written as 

XdV i + dV i - 1 =Z i dV n - 1 , i = 0,...,n-l, (9) 

where V-i :— 0. If we also let — det (X + fxl), the system ([5]) reduces 
to the system ([2j when we specify a metric on Q and consider the equiva- 
lent "vector version" A^VV^ =0 of ([5|). Restricting the attention to the case 
= dct (X + [il) is quite natural since, if we also assume that the coefficients 
Zq, . . . , Z n _i of are functionally independent, it is a necessary condition for 
the equation ([7]) to be satisfied. To see this, choose coordinates q 1 , . . . , q n as 
q 1 = Z^_\. Then the equation J7J), or equivalently XdZ t + dZ^_i = ZidZ n _i, 
reduces to Xj = q l 5™ — 5j~ , where S is the Kronecker delta symbol. In other 
words, —X must in these coordinates be the companion matrix (see (j!0|) ) of the 
polynomial Z^, and therefore it follows that Z^ = det (X + ptl). 

The *— multiplication reduces to Lundmarks multiplication of cofactor pair 
systems if we let X = J~ l J, where J and J are special conformal Killing tensors. 
In J4%, several other families of tensors X, that satisfy the equation (j4|) have been 
found. 

2.1 Matrix notation 

For the purpose of further study of systems of PDE's of the form ([5]), we 
shall introduce a new kind of matrix formulation for these systems and for 
the corresponding multiplication. The matrix formulation makes it possible 
to give an explicit formula for calculating powers of solutions, with respect to 
the *— multiplication. 
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The idea is to consider the column matrix V = [Vq, V\, . . . , V m -i] T instead 
of the polynomial = Vq + 1XV1 + • • • + /(" 1 " 1 V m -i, and to observe that 

V = V c e x = (V C° + ViC 1 + ■■■ + V m ^C m ^) e x 

with Vc = V^ = c where we have formally substituted the parameter \i with the 
companion matrix 



C := C[Z^} = 





1 

1 







o -z 

-Zi 



-Z n _ 2 

1 —Z n _i 



(10) 



of Z^, and where e\ — [1, 0, . . . , 0] T . The most important advantage of the 
matrix notation is that we can express the Euclidean algorithm of polynomial 
division in a more explicit way. For any polynomial = Pq + • • • + Pt/J? , 



the residue modulo Z^ can be written as R^ 



Rq + ■ ■ ■ + R m -\[i m \ where 



[Ro, R\, . . . , i? m _i] T = Pce-x- Thus, in the matrix notation, the *— multiplication 
of two solutions V and W of (O can be written as V * W — VcWc&x- In par- 
ticular, the matrix version of /x" is C a e\. The system ([5]) can also be expressed 
in terms of matrices as 

k m—1 

= E E ° l+3ei i 9iV j> • • • > dnV ^ Ai 

i=0 j=0 

fe /m-1 \ 

= E C M z2 ( "" dxVj,--- ,d n Vj] ^ 



i=0 
k 



3=0 



(11) 



i=0 



where we consider Ai as the matrix with elements (-Aj)jJ and V is the functional 
matrix 



V 



01 Vb 
diV m -i 



d n V 

driVm-l 



The equation pip is indeed independent of coordinates. The expression on the 
right-hand side of equation (jll[) isamX 1— matrix consisting of the 1— forms 
j{C r ) a jAi<lVj-i, a — 1, ...,m, where {C l ) a j is the element in the row a 
and the column j of the matrix C l . The * — multiplication theorem can then be 
expressed as: 
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Proposition 3 (Explicit criterion for existence of *— multiplication). Let V and 

W be solutions o/O = y^._ C l V ' A4, then V * W := VcWc&i is also a solution 
if and only if 



= J2 ClZ ' A * 



(12) 



i=0 



3 Explicit form of linear PDE's admitting 
*— multiplication 

The system (JS|) contains three parameters: n - the dimension of the manifold 
Q; m - the polynomial degree of the function Z^; k - the polynomial degree of 
the tensor A^. In this section we will discuss how the form of the system ([5]), or 
(jlip in matrix notation, and of the related *— multiplication depends on these 
numbers. We will also, for different choices of n,m,k, specify the structure of 
typical (generic) systems that allow *— multiplication. This is done by choosing 
the functions Zi as coordinates. 

When k — 0, the system (|S|) can be written as ^4odVi =0, i = 0, . . . , m — 1 
and the multiplication becomes trivial. When n = 1, the system ([S]) reduces to 
a quite simple system of ordinary differential equations. Therefore, we consider 
only cases where k > and n > 1. 



The *— product V * W of two solutions is a collection of m functions, each 
being a sum of functions of the form PijViWj, where Pij is a polynomial expres- 
sion of the variables Z^. The structure of the *— multiplication formula depends 
only on the parameter m, not on n or k. Since the degree of the polynomials 
Pij will not exceed m — 1, the multiplication will be more complex for higher 
values of the parameter m. For the simplest case m = 2, we have 

V * W = V c W c ei 

= (V W C Q + (V Q Wi + V1W0) C + V1W1C 2 ) e x 

V W ~Z V 1 W 1 
V Q Wx + VxW - ZxVxWx 

and for m = 3 the *— product V * W is given by 



V W - Z VxW 2 - Z Q V 2 Wx + Z Z 2 V 2 W 2 
V a Wi + ViWo - ZxVxWi - Z x V 2 Wx + {-Z + Z X Z 2 ) V 2 W 2 
V Q W 2 + V 1 W 1 + V 2 W - Z 2 VxW 2 - Z 2 V 2 W X + (-Zx + z$) V 2 W 2 



(13) 



3.1 Generic cases for different choices of (n,m,k) 

Our approach to find explicit forms of equations admitting *— multiplication, 
is to choose some generic coordinates in which the system (JSJ, equipped with 
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*— multiplication, takes a simple form. Since the functions Zi play a fundamen- 
tal role for the multiplication, we will assume that as many of these functions 
as possible are functionally independent and take them as local coordinates on 
Q. Since there are at most n different functionally independent functions, the 
relation between to (the number of functions Zi) and n will be crucial for spec- 
ifying each generic case. 

When to = n, we can choose as generic coordinates q 1 = Zq, . . . , q n = Z m -\ 
if the functions Z{ arc functionally independent. We consider the case when 
TO = n as our main case since the system (|5|) takes a simpler form than in the 
other cases. 



When to < n, the functions Zi are too few to form a complete set of co- 
ordinates. Instead we choose generic coordinates g 1 ,...,^™ such that q 1 = 
Zq, . . . , q m = Z m -i, without specifying the last n — m coordinates q m+1 , . . . , q n . 

When to > n, the functions Zi must be functionally dependent. For the 
generic case we assume that Zq, ...,Z n _x are functionally independent, and 
choose them as generic coordinates. 



We shall present below an explicit form of the system ([5]) in generic coor- 
dinates for the cases to = n, to < n, and to > n. For each case we will also 
consider simpler sub-cases according to the following schematic diagram (11411 
for the triples (n, to, k) 



(n,m<n,k) 



(n,m=n,k) (n,m>n,k) 



I 



I 



I 



(n,2,l), (n,3,2) (n,m=n,l 



I 

(3,2,1) 



I 

(2,2,1) 



(2,3,1) 



(14) 



For the sake of simplicity, we will consider the case when A m —i is non- 
singular. It is then no restriction to assume that A m -\ is the identity mapping. 

3.2 m = n 

A generic case with the simplest structure is obtained when to = n. According to 
the discussion above, (n,m,k) — (2,2,1) is the lowest value of the parameters 
for which the multiplication is non-trivial. It is also the best case to study 
in order to get a good understanding of the mechanism of the multiplication. 
We will investigate this case in detail, and after that some of the ideas will be 
generalized to the cases (n = to, to, k) and (n = to, to, k = 1). 
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3.2.1 (n, m, k) = (2, 2, 1) 



For this choice of parameters n, m and k, we have Z M = Z 
A + fiAi , = Vo + fiVi and ((3J) can be written as 

f A dV = ZoAidFi 

I AidVb = (ZiAi -A )dVi, 



\xZ Y + fi 2 



(15) 



or, as we have seen, in matrix notation as = V'Aq + CV A\. Thus, in local 
coordinates, (|15p constitutes a system of four partial differential equations for 
two unknown functions Voj Vi of two independent variables x, y. Since the num- 
ber of equations exceeds the number of dependent variables, this system will in 
general be overdetermined. If we assume that Ai is non-singular, we can instead 
of (|15[) consider the equivalent system 



= (A 2 - ZiA + Z Q I) dVi 
dV = {Z l I-A)dV ll 



(16) 



where A — AqA^ 1 , or in matrix notation = V'A + CV . We assume now 
that the functions x = Z$, y — Z\ are functionally independent and con- 
sider the system (116[) in the generic coordinates x,y. In these coordinates the 
relation (fP2|) . that guarantees a *— multiplication of the corresponding sys- 
tem (fT6| . reduces to A = —C. Thus, in generic coordinates, if we require 
existence of *— multiplication, the first equation of system (TlB|) is a conse- 
quence of the Cayley-Hamilton theorem so the system reduces to the system 
dVb = (yl + C) dVi, that has components: 



dVo 
dx 

m 

dy 



dV x dVi 
1 dx dy 
dVi 
dx 



Thus, the generic (2, 2, 1)— case constitutes in fact a determined system of two 
partial differential equations for two unknown functions of two independent 
variables. As we have seen, the *— product of two solutions can in this case be 
written as 



V * W = V c W c e 1 



VqWq - xViWi 
V W 1 + V 1 W -yV 1 W 1 
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The simplest non-trivial solutions obtained by taking powers of the trivial solu- 
tion (0, 1) are 



(o,i); 2 = 


( y 2 -x 

\ X 2 


(o,i)^ = 


(_»,_! 

V X X 


(0,1)2 = 


(1,0) 


(o,i): = 


(0,1) 


(o,i)2 = 


{-x,-y) 


(0,1)2 = 


(xy, -x + 



(17) 



y 



If the roots of the polynomial (or the eigenvalues of the companion matrix C) 
are functionally independent, we can instead define local coordinates through 
Zq = xy, Z\ = x + y. The condition ([7]) can then be expressed as Ao = DA\, 
where D = diag(a;, y). If we moreover assume that A\ is non-singular, the system 
([5]) reduces to DdVo = xydVi, or in components: 



( dV _ dVx 
dx ^ dx 
dVg _ dVx 

> dy dy ' 



This system has the general solution 

_ x(p(y) - yip{x) _ 4>{y) - ip{x) 

x — y x — y 

where <p and ip are arbitrary functions of one variable. 



(18) 



3.2.2 (n, m, k) — (m, to, k) 

We will now study the generic case of the more general situation when the 
only restriction for the parameters n, to, k is that n — to. We assume now 
that q 1 = Zo,...,q n = Z m -\ are functionally independent and constitute a 
complete set of coordinates. The condition (|12[) . which guarantees existence of 
multiplication for the system ([S]), now attains the simple form Ac = since Z' 
in these coordinates becomes the identity matrix. Hence, the system (J5J) admits 
*— multiplication if and only if it, in the matrix notation with the coordinates 
q l = Zi—i, can be written as 

k 

o = J2 ( civ> - v ' ci ) A * ( 19 ) 

i=l 

where the tensors A\,...,Aj i are arbitrary. 
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3.2.3 (n, m, k) — (to, to, 1) 

When k — 1, the equation P^|) becomes a remarkably simple equation = 
(CV — V'C)Ai in terms of the generic coordinates. Thus, if we also assume 
that A 1 is non-singular, we obtain in the generic (to, to, 1)— case, the equation 



CV' = V'C. 



(20) 



By calculating the residue of A^AV^ modulo Z^, we see that the equation ([20 
can also be written as 



-C&Vi + dVi 



.1 = q t+1 dV n . 



1, 2 = 0,1,. 



1, V-i := 



(21) 



However, the equation for which i = in (|2ip can be discarded since it is a 
consequence of the other equations and of the Cayley-Hamilton theorem for 
the companion matrix C. This is realized by adding to the equation in (|2ip for 
which i = the equation for which i = 1 multiplied with C, then adding the 
equation for which i = 2 multiplied with C 2 , and so on. In components, the 
equations (|2T|) can then be written as 



m 


dVi-i 






dq 2 

_m_ 


dq 1 
SVi-x 


= q l - 


dq 1 


dq 3 


8q 2 


= q l - 


dq 2 


dVi 






dV n -i 


dq n 

qn w_ 


g q n-l 


= q l - 


dq 11 - 1 


q dq n 


dq n 


= q l - 


dq n 



xdVi 

where i = 1, . . . ,n — 1. Thus, we see that the generic (to, to, 1)— case consists of 
to(to — 1) PDE's for to dependent variables, and is therefore an overdetermincd 
system when to > 1, that nevertheless has non-trivial solutions, e.g., Vi = Zi = 
q' +1 . If we let U = V * W, each entry U a of the 1— column matrix U is a sum 
of terms PijViWj, where Pij is a polynomial of degree at most to — 1 in the 
coordinates q , . . . , q n . We note also that for every solution V of equation (f2"0"l) , 
each term of the sum in equation (JT9J) vanishes. Thus, every solution V in the 
generic (to, to, 1)— case also solves the generic (to, to, A:)— case for arbitrary k. 
The lowest order *— powers of the trivial solution fj, (or [0, 1, 0, . . . , 0] T in matrix 



13 



notation) are given by: 

/4 = [0,1,0, 

A = [0,0,1, 



,of 
,of 



= W<r, -q 1 + ?V\ ■ • • , -r- 1 + (q m ) 2 ] T 

^+2 = [ q l {q m-l _ (9 m )2)j + q 2 {q ,n-l _ {q ,n )2 ^ _ 

3.3 m < n 

Suppose that Zq, . . . , Z m ^\ are functionally independent and consider local co- 
ordinates q , . . . ,q n such that q 1 — Zq, . . . , q rn = Z TO _i, without specifying the 
other coordinates q m+1 , . . . , q n . Now (fT2]) is equivalent to = J2i=o \p l I 0] 
where \C l | 0] denotes amy. n— matrix constructed by writing n — m zero- 
columns right to the matrix C '. This condition determines the first m rows of 
Aq uniquely, while the last rows as well as Ai, A 2 , . . . , A^ are arbitrary. 

One should note that when m < n, since Zi are the only functions except V* 
and Wi that appear in V * W, this product will not depend on certain coordi- 
nates unless V or W are themselves functions depending on these coordinates. 
Thus, a product of two trivial solutions will will never depend on these missing 
coordinates. 



3.3.1 (n, m, k) = (n, 2, 1) 



We will now choose the lowest possible values for m and k and let n be arbitrary. 
We note that the (3, 2, 1)— case, which is the simplest possible case for which 
m < n, is included. In generic coordinates, the relation (|12[) is equivalent to the 
following relation between the components of the tensors Aq and A\ : 



(A 
(A 



(A 
(A 



= C 



where C denotes the companion matrix of Z^. We see that, even though we 
restrict our attention to the case A\ — I, n(n — 2) components of A can still be 
chosen arbitrary. However, for some choices of these components, the system 
([5]) will not depend on some of the coordinates and can therefore be reduced to 
a lower dimensional problem with a smaller number of independent variables. 
We illustrate this phenomenon in the case (3, 2, 1). 
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3.3.2 (n, m, k) = (3, 2, 1) 



According to the discussion above, we can assume that 

"0x0 

A = -1 y 
a 6 c 



in generic coordinates (x, y, z) where a, b, c are arbitrary functions. By analyzing 
the corresponding system one can see that either all solutions Vb, V% are 
constant with respect to the variable z and the system then reduces to the 
generic (2,2, 1)— case, or otherwise we must have y 2 > 4x, b — a(c — y), and 
c = y ± (1/2) y/y 1 — Ax, and the system can then be written as AodVb = xdVx, 
or in components 



dV dV 8Vi 

- — + a(x,y,z)-^=x — 



dV dV 
x— +V-q^ +a{x,y,z) 





dVx 

' dy 



where a — a(x, y, z) is an arbitrary function. Thus, we see that the generic 
(3,2,1)— case involves an arbitrary function, which was not possible for the 
(m, m, l)-case. We note also that all solutions V of the generic (2, 2, l)-case 
also solve the generic (3, 2, 1)— case. Since the multiplication coincides with the 
multiplication in the (2, 2, 1)— case, we see that unless V or W depend on z, 
the product V * W will not depend on z. Especially, the solutions obtained by 
taking powers of the trivial solution fi are again given by (|17|) . 



3.3.3 (n, m, k) = (n, 3, 2) 

When we consider higher values of the number k, the corresponding systems 
become harder to analyze. One reason is that when we increase k by one, we 
add a new tensor A/, which means that we add n 2 new components. Another 
reason is that with higher values of k, we get higher order polynomials in Zj. 
Already for k = 2, such systems become quite hard to handle. In the general 
(n, 3, 2)— case for example, the system §5§ can be written as 

= A dV - Z A X &V 2 - Z A 2 dV 1 + Z Z 2 A 2 dV 2 

= Aod^i + AidVo - Z 1 A 1 dV 2 - Z^A^dVy + {Z X Z 2 - Z ) A 2 dV 2 

= A dV 2 + AidVi + A 2 dV - Z 2 A x dV 2 - Z 2 A 2 dV x + (Z 2 - Z x ) A 2 dV 2 . 

Even if we assume that A 2 = I and that q t+1 = Zi are functionally independent, 
we still have 2n 2 arbitrary functions in the picture (the components of Aq and 
Ai). When the condition p^]) . which in generic coordinates is equivalent to 
= Ei=i [C l I °] Ai, is satisfied we still have n(2n — m) arbitrary functions. 
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3.4 m > n 



When to > n, the generic case becomes more complicated than for m < n. 
Consider the generic case when q 1 = Zo,...,q n = Z n _\ are functionally in- 
dependent and take g 1 ,...,?™ as local coordinates. In these coordinates we 
have 



Z' = 



1 





diZ n 
diZ m -i 




1 



d 2 Z n 
d 2 Z m ~\ 






1 

d n Z n 

dn Zjri- 



Thus, the relation (fT2"| . which in the generic cases for to < n became a set of 
algebraic equations for the components of Ai , becomes now a complicated differ- 
ential relation between the components of Aj and the functions Z n , . . . , Z m _\. 

3.4.1 (n, m, k) = (2, 3, 1) 

We consider the simplest case for which to > n, i.e., when (n, m, k) = (2, 3, 1). 
The system ([5]) can then be written as 



= A dVb - Z AidV 2 
= AiAV a + A dVi - Z 1 A 1 dV 2 
k = A x 6Vi + {A a - Z 2 A 1 ) dV 2 . 



(22) 



The condition to > n, implies that the functions Zq, Z\, Z 2 are functionally 
dependent. For the sake of convenience, we assume in the generic case that the 
functions Zq, Z\ are functionally independent and that Z 2 = (^(Zo, Z\) for some 
function (f>. We also assume that A\ = /. Thus, in generic coordinates x = Zq, 
y = Zx, the condition (fTSj) is satisfied if and only if A := Aq is given by 



.4 



yd x cj> 



Xdy4> 

ydy<j) 



where <p satisfies the non-linear partial differential equations 



= x -M - 



dx 



dx ^ dx dy dy 



= 1 



Solutions of the equations 



d(j) dcj) 
dx dy 



dy 



dy' 



exist, for instance <f> = ay — a 2 x 



(|23D for any non-zero real constant a. In the generic case, the system 



(23) 



solves 



is 



1G 



equivalent to (note that the tensor A is non-singular when x =/= 0) 



= AdV - xdV 2 
= (xl - yA + cj)A 2 - A 3 ) dV 2 
= dVi + {A - <pl) dV 2 . 

The Cayley-Hamilton theorem together with the assumption that <f> satisfies the 
equations ([23]) implies that the expression xl — yA + <pA 2 — A 3, vanishes. Hence, 
the system (|^)) reduces to 




xdV 2 

{</>!- A) dV 2 . 



In the case when (f> = ay — a 2 x + a 1 , the system 
as the overdetermined system 



can explicitly be written 



2 9V 2 , 

ax— 1- (ay h x— — 



dx 
dx 



dV 2 
dx 



ay- 



dVo 
dy 

dy dy 

-u dy 2 . 2 ,.dV 2 
a 1 )^ + (a 2 y+D^ 



dv x dv 2 ! 2 ,sy 2 

-— ax— h (a —ax)——. 

ay 



dx 



dy 



For this system, the *— product of trivial solutions gives in general non-trivial 
solutions. For example, when V — [0,1, 0] T and W = [0,0,1] T (or simply 
— /i and = fJ, 2 using the fi— notation), their *— product becomes 



V*W = 



ax 



—x 

-y 

ay - 



3.5 Summary 

The results about the possible structures of the system of linear PDE's and 
of the corresponding relation ((7|) (that characterizes the existence of *— multi- 
plication), can for non-singular Ak be summarized as follows: 

1. to = n. Suppose that Zq, . . . , Z m -\ are functionally independent. In the 
generic coordinates q l = Zi-i, the relation ([7]) is satisfied if and only if 
Ac = AqC° + ■ ■ ■ AkC k = 0, and the system (|11[) can be written as 

k 

= Y J {C t V / -V'C^A,, 

i=l 
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where the tensors A\, . . . , Ak are arbitrary. In the generic case when k = 1, 
the system (fTTj) takes the remarkably simple form CV — V'C, for which 
every solution V is also a solution in each generic (n, to — n, ft)— case with 
arbitrary k > 0. 

2. to < ?i. Suppose that Zo, . . . , Z m -\ are functionally independent. In the 
generic coordinates q 1 = Z , . . . , q m = Z m -\i q m+1 , ...q n (q m+1 , ■ ■ ■ , q n 
not specified), the relation (J7]) is satisfied if and only if = Yli=o \p l I 0] ^i- 
It is characteristic for the case to < n that *— products of trivial solutions 
will remain constant in some variables. 

3. to > ft. This is the hardest case to analyze since the relation (O in the 
generic coordinates q 1 — Z$, . . . , q n = Z„_i becomes a differential relation 
for the functions Z n +±, . . . , Z m _i, while it leads to algebraic equations in 
the previous cases. 



4 Power series 

As mentioned above, we can take power series of the trivial solution /j, to build 
up more complicated solutions of ([5]). In the matrix notation, such power series 
have the form 



P = J2a r C r e 1 , 



r=0 

where a r are real constants. We will now investigate when these power series 
define new solutions of the system (J5]), i.e. when they are convergent and the 
summation and derivation commutes so that they define genuine solutions of 
the first order systems of PDE's ([5]). 

The companion matrix C can be factorized as C = TJT~ l , where the matrix 
J has the Jordan canonical form, i.e. J = diag(Ji, . . . , J s ), where 



X s 1 
A., 



1 

A. 



(24) 



and A s is an eigenvalue of C (note that the eigenvalues of C coincide with the 
roots of the polynomial Z^). Thus, the partial sums of the power series can be 
written as 



N / N \ 

p N = a rC r ei = T a r J r T -1 ei. 

r=0 \r=0 J 



(25) 
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Lemma 4. Let J s be the Jordan block defined by (I24[) . The entries of powers 
of J a can are given by 

(/,) A -" ' 

Proof. The proof follows immediately by induction over r. 

□ 

From the lemma, it is clear that all non-zero elements of the matrix X^-Lo a rJ r 
have the form J2r=o ar (b) K~ b i f° r eigenvalues X s of C. All power series 
ar (b) As~ b have the same radius of convergence as Y^Lo -^ '■ Hence, we 
conclude that at every point of Q where all eigenvalues of C belongs to the open 
set ] — R,R[ (where 1/R = limsup lani 1 /™), the power series ^2^L a r J r , and 
therefore also a rC r ei, will be convergent. The question that now remains, 

is what is required for a convergent power series X^lo a rC r e\ to be a solution 
to ([5]). The following theorem gives a sufficient condition: 

Theorem 5. Consider a power series P = X^Lo a rC r e\. Let D C Q be a 
domain in which each eigenvalue X s of the matrix C is real and X s G [-R + 
e,R— e] for some e > 0, where 1/R = limsup |a r | 1 / r . Assume also that the 
geometrical multiplicity of each eigenvalue is constant in D. Then the power 
series P defines a solution of ([3]) in D. 

Remark 3. The assumption about the geometrical multiplicity guarantees that 
the structure of the Jordan canonical form of C is preserved in D ( note that 
the case with simple eigenvalues is included). There is no indication that this 
assumption, or the assumption about eigenvalues being real, are necessary for 
the conclusion in theorem [5j but without them the proof would be technically 
more complicated. 

Proof. P is a solution of ([5|) if and only if 

k k / oo \ ' 

= C'P'A, = $^ C* I a r C r ei Ai. 

i=0 i=0 \r=0 ) 

Thus, since C r e\ is a solution for any r, we see that it is enough to prove that 
P' = liniTv Pjv, or in components, diP 3 — lim^r diP 3 N , i, j = 1, . . . ,n, where P 3 
and P J N denote the j'th element of P and Pjv, respectively. Since Pn converges 
to P in every point, we only have to prove that diP 3 N converges uniformly in D. 

According to (j25|) and lemma [H each P 3 N consists of a finite sum of terms 
of the form f(q) b r X r s , where X s is an eigenvalue of C, / is a real- valued 
function, and the power series X^ b r X T s has R as its radius of convergence. 
Thus, diP 3 N consists of a finite sum of terms of the form {dif(q))^2^ b r X r s + 
f(q)(diX s ) ^2^? rbrX^ 1 , which converge uniformly in D. □ 
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According to theorem [5j if the companion matrix C has simple eigenval- 
ues and J2 r a rt r i s a power series with infinite radius of convergence, then 
J2 r a r C r ei is a globally defined power series solution of Thus, for example, 
we can construct the power series solutions 




cos*C 



We end the discussion about power series with some examples. 

Example 3. To illustrate the mechanism of generating power series solutions, 
with respect to the * — multiplication, we return to the system (|15[) in the (2, 2, 1) — 
case. We assume that and satisfy the relation (J7J so that the system has 
* — multiplication, and we also assume that the roots Ai, A2, of are simple. 
The companion matrix can then be diagonalized as 



1 _ 


' A 2 Ai 




' Ai 







A2 Ai 




-1 -1 







A 2 _ 




-1 -1 



C = TDT- 



A power series ^2 r a r C r e\ defines a solution in any domain D C Q in which 
I Ai|, A2I < 1/limsup \a r \ x / r — e, and can be written as 



a r C r e 1 =T\S^ a r D r j T -1 ei 

r \ r / 



1 


A2 Ai 




' J2a r Xl 




" -1 -Ai " 




' 1 " 


Ai — A2 


-1 -1 




J2 a r^2 _ 




1 A 2 








1 



Ai — A2 
Thus, for example we have 

exp„C 



Ai ^2 ftj-A^ — A2 ^2 a r \\ 
a r Ai — Y2 a r A 2 



Ai — A2 



Ai exp (A 2 ) - A 2 exp (Ai} 
exp(Ai) -exp(A 2 ) 



In the case when the eigenvalues of C are constant ( which is the case for the 
Cauchy-Riemann equations), these power series will only provide trivial solu- 
tions. Hence, in order to obtain interesting solutions for this case, one has to 
build power series solutions from a non-trivial solution (for instance (x, y) in the 
Cauchy-Riemann case). The other extreme case is the generic situation when 
the eigenvalues are functionally independent and A\ is non-singular. Then, as 
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mentioned earlier, the general solution of (|18[) can be written as 

1 



V = 



\\ — A2 



Ai^(A 2 ) - A 2 V(Ai) 
V>(Ai)-0(A a ) 



(26) 



.ft is a remarkable property that every analytic solution in the generic case can 
be expressed by power series of the trivial solution (0,1). Namely, if V is a 
solution of the form (|26p where <p(t) = ^2 a r t r and ^(t) = ")2 b r t r are analytic, 
then V can be expressed in terms of the two simple solutions (Ai — A2) _1 (Ai, —1) 
and (Ai — A2) _1 (— A2, 1), and of power series of the trivial solution (0, 1): 



V = 



Ai — A2 



Ai 
-1 



arC 



Ai — A2 



-A 2 
1 



Hence, a significant part of the solution set of (|T8|) can be expressed by power 
series in trivial solutions. 

Remark 4. The solutions (Ai — A2)~ 1 (Ai,— 1) and (Ai — A2)~ 1 (— A2,l) have 
some remarkable properties. They are idempotent and their sum is the identity, 
i.e., 





Ai 


-1 


Ai 


-Aa 


' Ai - A 2 




-A 2 


1 


Ai 


-A 2 


Ai — A2 




Ai 


-1 


Ai 


-A 2 


' Ai — A2 



Ai 



-1 



Ai — A2 Ai — A2 

-A 2 1 

Ai — A2 Ai — A2 
-A 2 1 



Ai — A 2 Ai — Aj 



(1,0). 



Example 4. Let us instead consider the case where m = 2 and the eigenvalues 
of C coincide everywhere, Ai = A2 =: A, i.e., the functions Zq and Z\ are 
functionally dependent and related as Z\ = AZq. Then we can factorize C as 



1 


1 -A " 




' A 


" 




" 1 


A " 




1 




1 


A 







1 



and a general power series solution can be written as 
J2 a r C r e 1 =T[J2 a rJ r ) T- 1 ^ = 



J2 r (l-r)a r \ r 



For this example, the exponential power series produces the solution exp^C 
exp(A)(l - A, 1). 



Example 5. In J3]/ ; a matrix equation of the form 

V/ = MV S , 



(27) 



where M is a constant matrix, is studied in an open convex domain of a real or 
complex vector space. Here /(x) and <?(x) are scalar-valued functions defined 
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on the vector space, and V/ = [d Xl f, . . . , d Xn f] T , where x = \x\, . . . , x n ) T is 
a coordinate vector with respect some basis. We give below an example of an 
equation of the form (|27p that admits a * — multiplication, and we show that so- 
lutions can be represented by power series of simple solutions. 



Let V be a real vector space of dimension four, and assume that the matrix 
M in ([27]) has eigenvalues A, A, a±if3, where A, a, and (3 are real constants. By 
performing a linear change of variables x — > Ax., where A is a constant matrix, 
the equation (|27p transforms to V/ = A~ T MA T Vg. Thus, we can assume 
that M is in canonical real normal form. We assume that the eigenvalue A has 
geometrical multiplicity one, so that M is given by 



M = 



A 1 

A 

a (3 

-0 a 



(28) 



According to the general analytic solution of (|27|) can be decomposed as 
f = fi + f%, g = 91 + 32, where (/i, 31) and (/ 2 , g 2 ) are solutions of 



V/i 



A 1 
A 



V.gi 



id V/a = 



a (3 
—(3 a 



V32, 



(29) 



respectively. By changing to new dependent variables 
/i = /i-A3i, gi =01, h = f2-ctg 2 , 



h = f3g 2 , 



it becomes obvious that there is no restriction to assume that A = a = and 
(3=1. We note that the equation for / 2 and 32 then reduces to the Cauchy- 
Riemann equations. For this choice of M , it is trivial to obtain the general 
analytic solution of (|27|) : 



/1 = x 2 (f>'{xi) + ip(xi) 
31 = i>{xi) + c 



h = Re (F(x 3 + i Xi )) 
g 2 = lm(F(x 3 + ix 4 )) , 



where ip and <fi are arbitrary analytic functions of one real variable, and F is an 
arbitrary holomorphic function of one complex variable. We have already seen 
that the general solution of the Cauchy-Riemann equations can be represented as 
a power series of the simple solution X3 + 1x4 with respect to a * — multiplication 
(coinciding with multiplication of holomorphic functions). Also the first system 
in (|29|) admits a *— multiplication since it can be written as a system of the form 
(J5J) where 




-1 



fJ-I, V fl = fx + figi, 



and the relation |(7J) is trivially satisfied since is constant. It is remarkable 
that also for this system, the general analytic solution can written as a power 
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series of simple solutions. Namely, let the functions <j) and tp, from the gen- 
eral solution, be analytic with power series representations ip{ s ) — ^2 r a rS r and 
4>{s) — ^2 r b r s r , respectively. Then, the solution (/i,<?i) is given by 

fi + mx = ^2( a r + /AO * (xi + /*r 2 )*. 

r 

Thus, when M is given by (|28l) . we have shown that the general analytic solution 
of (1271) can be obtained by taking power series, with respect to * — multiplication, 
of simple solutions of the subsystems (|29[) . Also for general systems of type 
(|2"T|) . solutions can be built up by taking *— power series of simple non-trivial 
solutions. 



5 How to find systems with multiplication 

As we have seen, the problem of finding systems of the form ([5|) that allow a 
*— multiplication is equivalent to finding a (1, 1)— tensor and a function 
such that A^dZ^ = 0. In this general form, the problem is hard to handle, since 
in coordinates we may need to solve a system of complicated non-linear PDE's. 

One way to construct systems of PDE's allowing multiplication, is to choose 
the function Z^ first and treat it as fixed. The equation ([7]), or equivalently 
(fT2"|) , becomes then a system of linear algebraic equations for the components of 
Ai, which is easier to solve. It constitutes a system of mn equations for n 2 (k + l) 
unknown functions. The number n(k + 1) — m decides about how large family of 
solutions that can be found for each Z^. When n(k + 1) — m < 0, the number of 
unknown functions is less than the number of equations and we may not expect 
to find solutions for every choice of Z. On the other hand, for large values of 
n(k + 1) — m we get many families of solutions whenever the system is consis- 
tent. We illustrate this process of finding systems with a *— multiplication by 
an example. 



Example 6. Let Q be a manifold of dimension two with local coordinates x, y 
and suppose that Z M = x + y[i + xy^j 2 + fi 3 and k = 1. When Z^ is given, we 
can find all tensors A^ = Aq + [iA\ such that the corresponding system (JS|) is 
equipped with a *— multiplication, i.e, such that ([7]) is satisfied. 







Z'A 

1 



y 

o 
i 
o 



CZ'Ax 
(A))? 



-v 

-xy 





" 1 


" 







1 




. v 


X 



+ 



( A i)i 



(A)I 
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The solution of this system of linear equations is given by 

(x 2 y + l)f (x 2 y+l)g 



A = 



xy x 
•j 2 — 1 xy 



Ax, A x = 



x(l-y 2 )f x(l-y 2 )g 



where f = f(x,y) and g = g(x,y) are two arbitrary functions. The system ([5]) 
takes a form that is independent of the choice of f and g: 

tfl+^^W^-l)^ 00) 

o = (i + + .(i -v 2 )^ + - 3 d - y 2 )^+ 
(y 2 -i)(^ + i)^ 

Although there are two arbitrary functions f and <?, every choice give rise to the 
same system. Thus, for this particular choice of and parameters (n,m,k), 
(|30[) is the only system with * — multiplication. Since m = 3, the corresponding 
multiplication formula is given by (|13l) . 

In [3], the problem of finding systems with multiplication on Riemannian 
manifolds was studied for the special case when A^ = X + fil and = det A^. 
With those restrictions, the problem of finding a system equipped with a mul- 
tiplication reduces to finding a tensor X that satisfies the equation (U). The 
following families of solutions were found in [4 : 

1. X = J~ 1 J where J and J are arbitrary special conformal Killing tensors. 
In this case, the multiplication of cofactor pair systems [5] is reconstructed. 

2. Every X with a vanishing Nijenhuis torsion Nx = 0. This follows from the 
remarkable relation 2 (Xd(detX) - detXd(trX)) i = {N x )%C j k , where 
C = cof X, and that X and X + fil share the same torsion. This result 
holds also when no metric is specified on the manifold Q. 

3. If X is a non-singular solution, then X~ l is a solution as well. 

4. In [5], a method for constructing solutions X consisting of smaller blocks 
that satisfy ((4]) is presented. A similar result is valid for the more general 
equation 0: 

Theorem 6. Suppose that A^ and Z^ satisfies the relation (JJJ and let 
A^ be a (1,1) — tensor, and Z^ a function on a different manifold Q, sat- 
isfying the relation A^dZ^ = (mod Z^). Then the tensor A^ © A^ is a 
(1, 1) — tensor on the manifold Q x Q and it satisfies the relation 

(Afj, ® A^&iZ^Z^) = [mod Z^Z^). 
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Proof. 



(A M © A^d(Z^Z^) = (A M © 4) (dg(Z M Z M ) + d^Z^)) 
= A^Aq(Z ^) + A^&q{Z ^Z ^) = Z^A^AqZ^ + Z^A^dgZy, 
= ZpZ^a = (mod Z^Z^), 

where a is a certain 1— form, and Aq and dg denote the exterior differential 
operator on Q and Q, respectively. □ 

5. If the non-singular tensor X is diagonal in the coordinates q , ...,q n , 
then X satisfies the equation {4| if and only if X = diag(Xi, .. .,X S ), 
where X a — diag(0 a , ■ ■ ■ , 4>a) is a square diagonal matrix of size n a and 
(f>a = <t>a{q ni ^ . . , , h " a ) is an arbitrary (sufficiently regu- 

lar) function, depending only on the specified coordinates. 

As we have seen, there are several ways to construct system of partial differ- 
ential equations with multiplication of solutions. Nevertheless, for any given sys- 
tem of PDE's, it is in general hard to settle whether it admits a *— multiplication. 

6 Conclusions 

The *— multiplication constitutes a powerful method for generating, in a pure 
algebraic way, new solutions from known solutions of certain linear systems of 
PDE's. Especially by taking trivial solutions, we can construct large families 
of non-trivial solutions of systems for which non-trivial solutions are hard to 
obtain by other methods. By generalizing the ideas from [4], we have signifi- 
cantly extended the class of systems of PDE's admitting *— multiplication. By 
identifying which elements of the construction of systems of PDE's with multi- 
plication that are relevant, we have obtained much better understanding of the 
nature of the *— multiplication. Our insight into the mechanism of the multi- 
plication has been obtained due to the matrix formulation of the problem and 
of encoding the Euclidean algorithm as a matrix polynomial (section 12. 1|) , and 
due to the effective construction of power series solutions (section |4}. 

There are still many questions regarding the *— multiplication which it may 
be worth to study. Some examples of such questions are: 

1. Which solutions can be represented as power series of trivial solutions, and 
when do they constitute all solutions? What can we say about power series 
of non-trivial but simple solutions (compare with the Cauchy-Riemann 
equations where all holomorphic functions are represented by power series 
of a linear polynomial solution)? 
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2. In section[5l the problem of finding systems equipped with *— multiplication 
was treated. But the opposite problem is also interesting. Given a lin- 
ear system of PDE's, determine if it has a multiplicative structure on the 
space of solutions. 

3. Theorem Q] characterizes all systems that admit *— multiplication. We 
have also presented some typical (generic) systems with multiplication. 
However, in order to gain a better understanding of the systems that 
admit multiplication, it would be desirable to have some natural principle 
of classifying these systems. 

4. In example [5l we saw that a particular matrix equation of the form V/ = 
M\7g admits a *— multiplication and that the general analytic solution 
can be expressed by power series, with respect to the *— multiplication, 
of simple solutions. This is a much more general property, that can be 
generalized to every equation V/ = MV ' g where M is a constant matrix 
with either real or complex entries. The work on this problem is already 
in progress, and the results are being prepared for publication. 
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